Docking, Synthesis and Evaluation of the Antifungal Activity of Pyrimido [4,5-b]quinolins

In order to expand the application of Fe3O4@SiO2-SnCl4 in the synthesis of heterocyclic compounds, in this study, we wish to report the use of one-pot three component synthesis of pyrimido [4,5-b]quinolone derivatives (D1-D16) through reaction of 6-amino-2-(methylthio)pyrimidin-4 (3H)-one, dimedone, or 1,3-cyclohexadione and aldehydes in the presence of Fe3O4@SiO2-SnCl4 as an efficient eco-friendly catalyst under ultrasound irradiation. The final aim of this study is evaluation of antifungal activity of resulted products. Synthesis of pyrimido [4,5-b]quinolin derivatives were done via three components coupling reaction of aldehyde, dimedone or 1,3-cyclohexadione and 6-amino-2-(methylthio)pyrimidin-4 (3H)-one in the presence of Fe3O4@SiO2-SnCl4 under ultrasonic irradiation in water at 60 °C. The products structure were studied by FT-IRI, 1H NMR,II and 13C NMRII. All the compounds were screened for antimicrobial activity by broth microdilution methods as recommended by CLSIIII. Considering our results showed that compound (D13) had the most antifungal activity against C. dubliniensis, C. Albicans, C. Tropicalis, and C. Neoformance at concentrations ranging (MIC90) from 1-4 μg/mL. Compounds (D9), (D10), (D14), and (D15) had significant inhibitory activities against C. dubliniensis at concentrations ranging (MIC90) from 4-8 μg/mL, respectively. 5-(3,4-dihydroxyphenyl)-8,8-dimethyl-2-(methylthio)-5,8,9,10-tetrahydropyrimido[4,5-b]quinoline-4,6(3H,7H)-dione (D13) exhibited inhibitory and fungicidal activities against the tested yeasts. The specific binding mode or the binding orientation of more efficient compounds to CYP51 active site, have been also performed by molecular modeling investigations and showed that there is a good correlation with biological test.


Introduction
Pyrimidines are an important class of heterocyclic compounds present in various natural products and as scaffolds in numerous pharmaceutically important molecules and functional materials (1). Pyrimidoquinoline derivatives have been used in a number of biologically active compounds with anticancer (2), antimicrobial (3), antimalarial (4), and anti-inflammatory activities (5).
Few methods for the synthesis of pyrimido quinoline derivatives have been reported using MCRs in the presence of 1-n-butyl-3methylimidazolium bromide ([bmim]Br) (6), p-toluenesulfonic acid (p-TSA) ( NPs-cell) under ultrasonic irradiation (13) and under microwave irradiation (14). Owing to the wide applications and significance of pyrimido quinoline derivatives in organic synthesis, pharmacology, and industry, there is still the need to develop an efficient, green, mild, and environmentally benign protocol for the synthesis of these important compounds.
We have recently reported efficient and eco-friendly procedures for the preparation of xanthene and 1,4-dihydropyridine derivatives, using Fe 3 O 4 @SiO 2 -SnCl 4 (15,16) under ultrasound irradiation. In order to expand the application of Fe 3 O 4 @SiO 2 -SnCl 4 in the synthesis of heterocyclic compounds, in this study, we wish to report the use of one-pot three component synthesis of pyrimido [4,5-b] quinolone derivatives through reaction of 6-amino-2-(methylthio)pyrimidin-4(3H)one, dimedone, or 1,3-cyclohexadione and aldehydes in the presence of Fe 3 O 4 @SiO 2 -SnCl 4 as an efficient eco-friendly catalyst under ultrasound irradiation. Also, some of the synthesized compounds were evaluated for their antimicrobial activity. In addition, docking simulation was performed to show the binding mode of the most efficient compounds to CYP51 active site.

Material and methods
Chemicals and solvents were purchased from Merck and Aldrich companies. FT-IR spectra were recorded as KBr pellet on a Bruker, Equinox 55 spectrometer. 1 H NMR and 13 C NMR spectra were recorded on a 400 MHz Bruker DRX-400 in DMSO-d 6 as solvent and TMS as an internal standard. Melting points were obtained with a Buchi melting point B-540 B.V.CHI apparatus. Ultrasonic irradiations were done using Elmasonic S 40H ultrasonic cleaning.
General procedure for the synthesis of pyrimido [4,5-

b] quinolones
In a round-bottom flask, a mixture of 6-amino-2-(methylthio)pyrimidin-4(3H)one (1 mmol), aldehydes (1 mmol), dimedone or 1,3-cyclohexadione (1 mmol), and Fe 3 O 4 @SiO 2 -SnCl 4 (0.03 g) in distilled H 2 O was sonicated in an ultrasonic cleaning unit (Elmasonic S 40H) at 60 °C for the stipulated time mentioned in table 2. The progress of the reaction was monitored by TLC (EtOAc:petroleum ether 7:3). After completion of the reaction, the catalyst was separated by an external magnet and reused for the next experiment. The reaction mixture was cooled to room temperature and then poured in to cold water. The solid product was filtered and washed with boiling water and recrystallized from ethanol to give the pure product in excellent yield.

Determination of antifungal activities Microorganisms
The antifungal activities of the synthetic compounds against some American Type Culture Collection (ATCC) strains of fungi,

Determination of minimum inhibitory concentration
MICs were determined using the broth microdilution method recommended by the CLSI with some modifications. Briefly, for determination of antimicrobial activities against fungi, serial dilutions of the synthetic compounds (1-1024 μg/mL) were prepared in 96-well microtiter plates using RPMI-1640 media (Sigma, St. Louis, MO, USA) buffered with MOPS (Sigma). Stock inoculums were prepared by suspending three colonies of the examined yeast in 5 mL sterile 0.85% NaCl, and adjusting the turbidity of the inoculums to 0.5 McFarland standards at 530 nm wavelengths (this yields stock suspension of 1-5 × 106 cells/mL). For moulds (Aspergillus spp. and Dermatophytes), conidia were recovered from the 7-day old cultures grown on potato dextrose agar by a wetting loop with tween-20. The collected conidia were transferred in sterile saline and their turbidity was adjusted to OD = 0.09-0.11 that yields 0.4-5 × 106 conidia/mL. Working suspension was prepared by making a 1/50 and 1/1000 dilution with RPMI of the stock suspension for moulds and yeasts, respectively. Working inoculums (0.1 mL) were added to the microtiter plates, which were incubated in a humid atmosphere at 30 ºC for 24-48 h. Uninoculated medium (200 μL) was included as a sterility control (blank). In addition, growth controls (medium with inoculums but without antibiotics or the synthetic compounds) were also included. The growth in each well was compared with that of the growth in the control well.
Docking procedure An in house batch script (DOCK-FACE) for automatic running of Auto Dock 4.2 was used to carrying out the docking simulations in a parallel mode, using all system resources as described before. To prepare the receptor structure, the complex of Mycobacterium tuberculosis-CYP51 enzyme with Fluconazole (PDB ID: 1EA1) was acquired from Protein Data Bank (PDB data base; http://www.rcsb. org) and water molecules and co-crystal ligand were removed from the structure. The PDB were then checked for missing atom types with MODELLER 9.17.
The ligand structures were made by using HyperChem software package (Version 7, Hypercube Inc). For geometry optimization, Molecular Mechanic (MM + ), followed by semi empirical AM1 method was performed. The prepared ligands were given to 100 independent genetic algorithm (GA) runs. 150 population size, a maximum number of 2,500,000 energy evaluations and 27,000 maximum generations were used for Lamarckian GA method. The grid points of 40, 40, and 40 in x-, y-, and z directions were used. A grid spacing of 0.375 A° was built centered on Hem group in the catalytic site of the receptor. Number of points in x, y, and z was -18, -3, and 67, respectively. All visualization of protein ligand interaction was evaluated using VMD software.

Results and Discussion
Synthesis Fe 3 O 4 @SiO 2 -SnCl 4 as an efficient nano magnetic catalyst is a readily available catalyst. This catalyst is removed from the reaction medium by an external magnet. In our ongoing research for new synthetic methods for the development of efficient and environmentally friendly protocols  Table 1). The products were characterized by their melting points using FT-IR, 1 H, and 13 C NMR spectroscopy.
The reusability of Fe 3 O 4 @SiO 2 -SnCl 4 was examined for four times without any considerable decrease in its efficiency.
The mechanism of the reaction and formation of compound 4a can be explained by the condensation, addition, cyclization, and dehydration reactions. The plausible mechanism for the synthesis of pyrimido [4,5b]quinolones in the presence of Fe 3 O 4 @SiO 2 -SnCl 4 , which can act as Lewis acid catalyst is depicted in scheme 2.
The carbonyl oxygen of aldehyde coordinates with the Lewis acid moiety increasing the electrophilicity of the carbonyl carbon and thereby making it possible to carry out the reaction in short time. In a plausible mechanism, it is assumed that the reaction may proceed initially through the Knoevenagel condensation between aryl aldehydes and dimedone or 1,3-cyclohexadione to form intermediate (I also act as a mild base for the deprotonation of an acidic proton of Intermediate (II). Then, intermediate (III) converts to (IV) via cyclization. Finally, the desired product (V) is obtained after dehydration of (IV) (Scheme 2).

Antifungal activities of the synthetic compounds
In this study, compounds (D1-D16) evaluated against fungi ( Table 2). None of the compounds have any effect on bacteria.
Considering our results showed that compound (D13) had the most antifungal activity against C. dubliniensis, C. Albicans, C. Tropicalis, and C. Neoformance at concentrations ranging (MIC90) from 1-4 μg/mL. Compounds (D9), (D10), (D14), and (D15) had significant inhibitory activities against C. dubliniensis at concentrations ranging (MIC90) from 4-8 μg/mL, respectively. In comparison of the antifungal activities of the synthetic compounds based on variation of substitutions on 2,3, and 4-position of phenyl ring, we found that the compound D13 with OH residue in meta and para positions of phenyl ring exhibited a better antifungal activities against the tested fungi than the other compounds. Compounds (D9) and (D10), had significant inhibitory activities against C. Dubliniensis, respectively. we found that the compound (D9) with Br residue in para position of phenyl ring exhibited a better antifungal activity against the C. Dubliniensis than the compound (D10) with CH 3 residue in para position of phenyl ring.

Scheme 2.
A proposed mechanism for preparation of pyrimido [4,5-b]quinolones in the presence of Fe 3 O 4 @SiO 2 -SnCl 4 Table 2. Minimum inhibitory and fungicidal concentrations of the synthetic compounds (µg/mL) against the examined fungi.  Table 2. Minimum inhibitory and fungicidal concentrations of the synthetic compounds (µg/mL) against the examined fungi Docking study of active compounds Docking is a method, frequently used to find the binding orientation of small molecule drug candidates to their protein targets in order to predict and interpret the affinity and activity of the small molecule. We performed molecular docking studies on D9, D10, D13, and Fluconazole to find, compare and validate their binding site, binding modes, and their best direction according to their binding energy (Figures 1).
Once the docking procedure was completed, to find the type of interactions, the protein-ligand complex was studied. All the docking protocols were done on validated structures, with RMSD values below 2 Å. The conformation with the lowest ones was considered as the best docking result. Docking binding energies of these active compounds were summarized in Table  3. Our results indicated that overall there is a good correlation between experimental pIC50 and docking binding energy ( Figure 1). As displayed in Table 3, all investigated complexes showed better docking binding energies than the co-crystal ligands (fluconazole). Consistent with our biological results, D13 was observed to have the best docking binding energies with the receptor through polar interaction with phe 255, His 259, and Thr 260 which responds to the key ineractions of pharmacophoric elements in the ligands (Figures 1). The docking model also indicated the specific alignment of fused three ring system to the pocket of heme iron of CYP51 in each ligand with similar positioning of triazole ring of fluconazole. This orientation with respect to the hydrogen bonding and hydrophobic interactions may be in favor of antifungal activity.

Conclusion
In conclusion, a convenient and efficient method for the synthesis of pyrimido [4,5-b] quinolin derivatives in the presence of Fe 3 O 4 @SiO 2 -SnCl 4 as a novel heterogeneous solid acid catalyst has been developed through ultrasound irradiation. This reaction under ultrasound irradiation not only gave excellent yield of products with lesser reaction time but it also required mild reaction conditions and exhibited operational simplicity. In addition, low cost, availability, recyclability, low toxicity of the catalyst and the use of water as a green reaction medium make this methodology a valid contribution to the existing processes in the field of one-pot multicomponent reaction.